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I. INTRODUCTION 

The momentum correlations of two or more particles at 
small relative momenta in their center-of-mass (cm.) system 
are widely used to study space-time characteristics of the pro- 
duction processes on a level of fm = 10^'^ m, so serving as a 
correlation femtoscopy tool. Particularly, for non-interacting 
identical particles, like photons or, to some extent, pions, these 
correlations result from the interference of the production am- 
plitudes due to the symmetrization requirement of quantum 
statistics (QS) fliSSl. 

The momentum QS correlations were first observed as an 
enhanced production of the pairs of identical pions with small 
opening angles (GGLP effect 01]). Later on, Kopylov and 
Podgoretsky settled the basics of correlation femtoscopy in 
more than 20 papers (see a review |5]) and developed it as a 
practical tool; particularly, they suggested to study the inter- 
ference effect in terms of the correlation function, proposed 
the mixing techniques to construct the uncorrected reference 
sample and clarified the role of the space-time characteristics 
of particle production in various physical situations. 

There exists |0,[llH|5lSI3.[ll an analogy of the momen- 
tum QS correlations of photons with the space-time correla- 
tions of the intensities of classical electromagnetic fields used 
in astronomy to measure the angular radii of stellar objects 
based on the superposition principle (HBT effect) ifioll . This 
analogy is sometimes misunderstood and the momentum cor- 
relations are mixed up with the space-time (HBT) correlations 
in spite that their orthogonal character and thus the absence 
of the former in astronomy measurements due to extremely 
large space-time extent of stellar objects (and vice versa) was 
already pointed out in early paper 1^] (see also |9]). 

Note that though the space-time (HBT) correlations are ab- 
sent in subatomic measurements, they can still be used in lab- 
oratory as an intensity-correlation spectroscopy tool allow- 
ing one to measure the spectral line shape and width (see 
iflll [T2I1 and references therein). In fact, the QS space-time 
correlations give the information about shape and width of 
the three-momentum distribution (including the angular one) 
of the quanta coming from a distant source and can be consid- 
ered as a spectroscopy tool in this more general sense. 

The momentum correlations of particles emitted at nuclear 
distances are also influenced by the effect of final state interac- 
tion (FSI) ill|Il[Tl[M[il- Thus the effect of the Coulomb 
interaction dominates the correlations of charged particles at 
the relative momenta 2k* in the two-particle rest frame smaller 



than the inverse Bohr radius |fl|^^ of the two-particle system, 
respectively suppressing or enhancing the production of parti- 
cles with like or unlike charges. 

Though the FSI effect complicates the correlation analy- 
sis, it is an important source of information allowing for the 
coalescence femtoscopy (see, e.g., 1 18, 19, 20, 21]), the corre- 
lation femtoscopy with unlike particles l,15> .16. ,17.1 including 
the access to the relative space-time asymmetries in particle 
production illlllMlHlllll^llillilMIIllIlIll 
and a study of strong interactions between specific particles 

In fact, the femtoscopic Coulomb correlations were ob- 
served more than 70 years ago when a sensitivity of the dif- 
ferential beta-decay rate to the nucleus charge and radius has 
been established [34]. 

The paper is organized now as follows. In sections II-IV 
we give the basics of formalism of femtoscopic and spectro- 
scopic correlations. The similarity and difference of femto- 
scopic correlations in multiparticle production and beta-decay 
is discussed in section V. We conclude in section VI. For re- 
cent results from femtoscopy of multiparticle production pro- 
cesses, one can inspect reviews 11331 ISSllSalSTI] . 



n. FEMTOSCOPIC QS CORRELATIONS 
A. Formalism 

The correlation function 2^(pi,/?2) of two particles with 
four- momenta p\ and p2 is usually defined as the ratio of 
the measured momentum distribution of the two particles to 
the reference one obtained by mixing particles from differ- 
ent events of a given class, normalized to unity at sufficiently 
large relative momenta. 

For identical pions or kaons, the effect of the strong FSI 
is usually small and the effect of the Coulomb FSI can be in 
first approximation simply corrected for (see |38] and refer- 
ences therein). The corrected correlation function is then de- 
termined by the QS symmetrization only. The space-time in- 
formation contained in the momentum QS correlations can be 
extracted based on the Kopylov-Podgoretsky (KP) formalism 
of independent one-particle classical emitters that are charac- 
terized by their production (excitation) four-coordinates, four- 
velocities, lifetimes and decay amplitudes. The KP formalism 
is valid on the following assumptions; 

(i) The multiplicity of produced particles is assumed suffi- 
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ciently large to neglect the influence of conservation laws. 

(ii) An independent or incoherent particle emission is as- 
sumed, i.e. the coherent contribution of multiparticle emitters 
is neglected (this contribution can be eventually taken into ac- 
count with the help of the suppression parameter X). 

(iii) The mean freeze-out phase space density is assumed 
sufficiently small so that the correlation function of two parti- 
cles emitted with a small relative momentum Q = 2k* in their 
cm. system is mainly determined by their mutual correlation. 

(iv) The momentum dependence of the one-particle emis- 
sion probabilities is assumed inessential when varying the par- 
ticle four- momenta pi and p2 by the amount characteristic for 
the correlation effect. This so-called smoothness assumption 
requires the components of the mean space-time distance be- 
tween particle emitters much larger than those of the space- 
time extent of the emitters. 

The probability amplitude to observe a particle with the 
four-coordinate x from an emitter A decaying at the four- 
coordinate XA can depend on x through the relative four- 
coordinate x—xa only and so can be written in the form: 



(x|\|/a) = (27t) *J c/"^K MA(K)exp[/K(x — jca)]. 



(1) 



We assume here that after the production the emitter moves 
along a classical trajectory and decays according to the expo- 
nential decay law. Such a classical treatment of the decay is 
often applied also to resonances (see, e.g., lf39l l40ll4lil ). It is 
justified for a heavy emitter with the energy spectrum of the 
decay particle substantially wider than the emitter width. One 
can avoid the classical treatment of the emitter decay by con- 
sidering xa in Eq. ([T]i as the emitter production (excitation) 
four-coordinate and adding, in the case of a negligible emitter 
space-time extent, the theta-function 9(f — Ja)- 

The probability amplitude to observe the particle with the 
four-momentum p is 



d^x {p\x){x\^a) = MA(/?)exp(-i>jt:A), (2) 



where {p\x) — exp(— The probability amplitude to ob- 
serve identical spin-0 bosons with four-momenta p\ and p2 
emitted by emitters A and B should be symmetrized in accor- 
dance with the requirement of QS: 



'AB 



\P1.P2) = [{piWa){p2Wb) + {p2Wa){p\W^) 



(3) 

The corresponding QS correlation function is given by the 
properly normalized square of this amplitude averaged over 
all characteristics of the emitters: 

1i.{P\,P2) = l 

, 3^ Lab ma (p 1 ) Mb (/52 ) m1(p2 )mb (f 1 ) exp (- /g'Ax) 



1 



Lab|ma(/?i)mb(;'2)P 
(cos(<7Ax)), 



(4) 



where p\ — P2, = xa —xb and the last equahty follows 
from the smoothness assumption (iv). 

It should be noted that the last equality in Eq. dUi is usually 
used to calculate correlation functions within classical trans- 
port models identifying the emitter four-coordinates as those 



of the decay four-coordinates of the primary emitters includ- 
ing resonances or those of the last collisions of the emitted 
particles. Concerning the accuracy of the classical approach 
to the emitter decay, we note that, for example in the case of a 
p-meson and a pion emitted from a small space-time region, it 
is a valid approximation at Q < 0.1 GeV/c but overestimates 
the tail of the corresponding two-pion correlation function by 
about 15 percent (see figure 2 in [39]). This overestimation is 
however not important when the space-time separation of the 
emitters is larger than 2 fm (as in heavy ion collisions) and so 
the interference effect rapidly vanishes at 2 > 0.1 GeV/c. 

It is instructive to introduce the emission function G{x,p) 
(similar to Wigner function) as a partial Fourier transform of 
the space-time density matrix: 

d'^E exp(-/pe)^(i + -e|\(/A)(VA|-«-2e>- (5) 

Since the single-particle emission probability is given by the 
integral 



d'^x G{x,p) 



El(/'l¥A> 

A 



(6) 



the emission function, though not positively defined, can be 
usually interpreted as a probability density to find a particle 
with four-momentum p and an average four-coordinate x = 
j{x+x'). For the QS correlation function of two identical 
spin-0 bosons one has: 



^{PUP2) = 1 



1 



/ d^x\d'^X2 G{xi , p)G{x2, p) cos{qAx) 

J d'^x\d^X2 G{xi,pi)G{x2,p2) 
(cos(^Ai)), (7) 

where Ax ~ xi — X2 and p = j{p\+ P2) = \P- Similar to Eq. 
dUl, the last equality follows from the smoothness assumption 
allowing one to identify the average four-coordinate x of the 
emitted particle with the decay four-coordinate of its emitter, 
i.e. neglect the space-time extent of the emitter 

Generally, for two identical bosons or fermions with the 
total spin S, the "H-" sign in Eq. ((Sj should be substituted 
by (^1)'^ or, equivalently, the two-particle ampUtude should 
be symmetrized (antisymmetrized) only for identical bosons 
(fermions) emitted with the same spin projections. As a result, 
in the case of initially unpolarized spin-y particles, the 
sign in Eq. (|4| or O for the QS correlation function should 
be substituted by (— l)^-'/ny, where nj is the number of pos- 
sible spin projections: = 2j + 1 or, for massless particles, 
nj = 2. 

A special comment requires the correlation function of 
massless particles with the spin j > 1/2 when the helicities 
in-between the extreme values ± j are forbidden. Thus, for 
photons, due to the transversality of the electromagnetic field 
(forbidden zero helicity), the spin factor multiplying the cor- 
relation term (cos(^Af)) actually depends on the angle 9 be- 
tween the photon three-momenta 143,143(1 . Particularly, in the 
model of randomly oriented one-photon dipole or quadrupole 
emitters, it equals ^(1 +cos^9) or ^(1 — 3cos^9 + 4cos'*9), 
respectively I43H . However, in the case of practical interest. 
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when the photon wavelength is essentially smaller than the 
size of the photon emission region, the angle between pho- 
tons contributing to the interference effect is very small and 
the spin factor reduces to the value 1 /tij = j irrespective of 
the multipole order of the emitter 

A comment requires also the correlation function of neutral 
kaons. Two neutral kaons with four-momenta pi and p2 are 
originally produced as pairs K'^{pi)K^{p2), K^{pi)K^{p2), 
K^\pi)K^{p2) and K^{pi)K^{p2). The coiTelation pattern 
now depends on the way the neutral kaons are detected. In 
principle, one can detect and K^, e.g., by charge exchange 
reactions and K . In this case, only the first 

two production channels would give the interference effect 
similar to that for identical pions. Usually, however, neutral 
kaons are detected through their two-pion and three-pion de- 
cays as so-called short-lived (/T^) and long-lived (Kf) states, 
respectively. Neglecting the small effect of Cf -violation, 
these states correspond to the eigen states with CP = ±1: 

\kI) = {\K') + |^°))/V2, \Kl) = {\K') - |^0))/V2. (8) 

Therefore, all four production channels of the pairs of 
K^- and ^"-mesons contribute to the production of the 
pairs Klipi)Klip2), Kl{pi)Kf^{p2), K^sipi)Kl{p2) and 
Kl{pi)K^^{p2). It is clear that the channels K^{pi)K^{p2) 
and K^{p\)K'-'{p2) yield the constructive interference pattern 
for any combination of the detected pairs of K'I- and K^- 
mesons. It is interesting that also the channels K^{pi)K^\p2) 
and K^\pi)K^{p2) now yield the interference effect which can 
be both constructive and destructive. Indeed, since there are 
two indistinguishable amplitudes contributing to final state, 
the symmetrized amplitude has the form [i441 : 

7xr(^f(/^i),^°(P2)) 

= {Kf\K'){K'j\K'){p,\y^^{K')){p2\^^^{K'^)) 
+ {k'^\K^){k'^\k'){P2Wa{k'')){p,Wb{K^)), (9) 

where the sign "H-" corresponds to Bose statistics of neutral 
kaons. One may see that the corresponding interference pat- 
tern is constructive for K^K^- and K^K^-pairs while it is de- 
structive for KgK^-pairs. 



B. Simple parameterizations 

A characteristic feature of the QS correlation function of 
two identical bosons (fermions) is the presence of the inter- 
ference maximum (minimum) at small components of the rel- 
ative four-momentum q with the width reflecting the inverse 
space-time extent of the effective production region. For ex- 
ample, assuming that for a fraction X of pion pairs, the pi- 
ons are emitted independently according to one-particle am- 
plitudes of a Gaussian form characterized by the space-time 
dispersions and Tq while, for the remaining fraction (1 — A,) 
related to very long-lived emitters (r|, r|', K^, A, . . . ), the rel- 
ative distances r* between the emitters in the pair cm. system 



are extremely large, one has 

^iPUPi) = l+A.exp(-r^q2-Tg^g) 

= 1 +Xexp{-rWT-{>-l + v\l)ql) , (10) 

where qj and qi are the transverse and longitudinal compo- 
nents of the three-momentum difference q with respect to the 
direction of the pair velocity \ — P/Pq. One may see that, due 
to the on-shell constraint [3] = vq = vqi (following from 
the equality qP — 0), strongly correlating the energy differ- 
ence ^0 with the longitudinal momentum difference qi, the 
correlation function at vZq > fQ substantially depends on the 
direction of the vector q even in the case of a spherically sym- 
metric spatial form of the production region. 

Note that the on-shell constraint makes the ^-dependence 
of the correlation function essentially three-dimensional and 
thus makes impossible the unique Fourier reconstruction of 
the space-time characteristics of the emission process. Par- 
ticularly, in pair cm. system, q — {0,2k*}, Ax = {t*,r*}) 
and the scalar product qAx = — 2k*r* is independent of the 
time difference t*. However, within realistic models, the di- 
rectional and velocity dependence of the correlation function 
can be used to determine both the duration of the emission and 
the form of the emission region [3], as well as - to reveal the 
details of the production dynamics (such as collective flows; 
see, e.g., ||45l|4^ and the reviews 10, HH]). For this, the cor- 
relation functions can be analyzed in terms of the out (x), side 
(y) and longitudinal (z) cornponents of the relative momen- 
tum vector q = {^.i,^vj?z} li49l Isoll : the out and side denote 
the transverse, with respect to the reaction axis, components 
of the vector q, the out direction is parallel to the transverse 
component of the pair three-momentum. The corresponding 
correlation widths are usually parameterized in terms of the 
Gaussian correlation radii Ri, 

ilipuPl) = 1 +Xexp{-Rlql-Rlql-R^^q^^-Rlq,q,) 

(11) 

and their dependence on pair rapidity and transverse momen- 
tum is studied. The form of Eq. (fTTT i assumes azimuthal sym- 
metry of the production process |47, 49]. Generally, e.g., in 
case of the correlation analysis with respect to the reaction 
plane, all three cross terms q/qj contribute ifsill . 

It is well known that particle correlations at high energies 
usually measure only a small part of the space-time emission 
volume, being only slightly sensitive to its increase related to 
the fast longitudinal motion of particle emitters. In fact, due 
to limited emitter decay momenta of few hundred MeV/c, the 
correlated particles with nearby velocities are emitted by al- 
most comoving emitters and so - at nearby space-time points. 
The dynamical examples are resonances, colour strings or hy- 
drodynamic expansion. To substantially eliminate the effect 
of the longitudinal motion, the correlations can be analyzed 
in terms of the invariant variable Q = {—q^Y^^ = 2k* and the 
components of the momentum difference in pair cm. sys- 
tem (q* = Q = 2k*) or in the longitudinally comoving system 
(LCMS) 115211 ■ In LCMS each pair is emitted transverse to the 
reaction axis so that the relative momentum q coincides with 
q* except for the component q^ — Jtll, where y, is the LCMS 
Lorentz factor of the pair. 
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Particularly, in the case of one-dimensional boost invariant 
expansion, the longitudinal correlation radius in the LCMS 
reads 1461 ~ [T /nitYl^i, where T is the freeze-out tem- 
perature, T is the proper freeze-out time and is the trans- 
verse particle mass. In this model, the side radius measures 
the transverse radius of the system while, similar to Eq. ( fTOl i. 
the square of the out radius gets an additional contribution 
{pt jnitf' IS.'^ due to the finite emission duration Ax. The ad- 
ditional transverse expansion leads to a slight modification of 
the /7(-dependence of the longitudinal radius and - to a no- 
ticeable decrease of the side radius and the spatial part of the 
out radius with p,. Since the freeze-out temperature and the 
transverse flow determine also the shapes of the -spectra, 
the simultaneous analysis of correlations and single particle 
spectra for various particle species allows one to disentangle 
all the freeze-out characteristics i47ll . 

III. SPECTROSCOPIC QS CORRELATIONS 

To help in understanding the analogy and difference of the 
QS space-time (spectroscopic) and momentum (femtoscopic) 
correlations, we briefly present here the formalism of QS cor- 
relation spectroscopy within the KP model of independent 
single-particle emitters. In the spectroscopic correlation mea- 
surements the particles are supposed to be emitted by a dis- 
tant object with large space-time dimensions and detected by 
two detectors at space-time points x\ and x^. It is assumed 
that the distance between the detectors is much smaller than 
the size of the emitting object and that this size is negligible 
compared with the distance between the object and detectors. 
Then the four-momentum of a photon emitted by the emitter 
A and detected by any of the two detectors can be written as 
PA = ft)A{l,/3A}, where /3a is the unit vector in the direction 
from the emitter A to the detectors. The four-dimensional in- 
tegral in the single-photon probability amplitude in Eq. ^ 
then reduces to the one-dimensional one: 

(x|\|/a) = (27l)"'y c/COa MA(C0A)exp[//:>A(x-Ji:A)]9(f-fA), 

(12) 

where maCwa) °= (wa — coqa — 'I^a) ^ i-C- the emitter decay 
is treated quantum-mechanically and parameterized by the en- 
ergy (Boa and width Fa of the emission hne. In accordance 
with the comment after Eq. ([T]i, now denotes the excitation 
four-coordinate of the emitter and the condition f > fA is intro- 
duced by the theta-function. Since the time fA is distributed in 
a very wide interval, the sum LaSxp[— /(coa — W^)fA] yields 
the delta-function 5((Ma ~ K)^) so that the single-photon prob- 
ability is merely proportional to the integral of the spectral 
function: 

EK-^IVa>P-E /^«A |«a(coa)|'. (13) 

A A 

The probability amplitude of two photons with the same and 
complete polarization should be symmetrized similar to Eq. 
©: 

'r^"(xi,X2) = [(jci|¥a)(jc2|v|/b) + (jc2|¥a)(jci|¥b)]/V2. 

(14) 



For photons with polarization f , the symmetrized amplitude 
(fl4l l describes only the fraction ^(l of the photon pairs. 
As a result, the correlation function 'S^(x\,x-i), defined as a 
number of two-photon counts normalized to unity at a large 
space-time separation of the detected photons, becomes 

1 +p2 

■s{,{xx,x2) = \^ — —■ 

LAB/^K'A^(0B|MA(C0A)P|MB((0B)pCOs(gABA)i:) 

Lab / -^Wac/wb I MA (coa ) n MB (cob ) P 
1 +p2 

= 1 + — ^ — (cos(^abAx)), (15) 

where q'ab =PK- Pb, Ax = xi- X2. 

It should be noted that the HBT technique is not based on 
counting separate quanta, it overcomes this difficult problem 
by the measurement of the product of fluctuating parts of the 
electric currents from the two detectors (the low-frequency 
part is filtered out) integrated and read in time intervals of the 
order of minutes. To get rid of the uncertainties in the detector 
gains the measured quantity is the ratio of the product mean 
to root-mean-square deviation. It can be shown that the cor- 
relation coefficient defined as this ratio normalized to unity at 
zero distance between the detectors, equals to 

(cos(qABAx)>»[2yi(p)/p]l (16) 

The approximate equality is valid for a uniformly radiating 
disk with the normal directed to the detectors (or, for a spheri- 
cal surface radiating according to Lambert law) emitting light 
of a small band width; the argument of Bessel function Ji is 
p = (odd, where CO is the mean angular frequency (mean en- 
ergy of the detected photons), d is the distance between the 
detectors perpendicular to the direction to the distant emitting 
object and is the object angular radius. Measuring the cor- 
relation coefficient as a function of the distance d, one thus 
determines the transverse spread Co9 of the wave vectors of 
the detected light (the spread of the transverse photon mo- 
menta). Obviously, the HBT correlation effect is insensitive 
to the actual space-time extent of the source. At most, one 
can determine the source angular radius 9 performing the ad- 
ditional spectral measurements of the mean angular frequency 

CO. 

It is interesting to note that Eq. ( fT6] l follows also from the 
superposition principle applied to classical electromagnetic 
fields and so the HBT intensity correlation effect would sur- 
vive even in the case of vanishing Planck constant when the 
QS correlations become unobservable. 

Comparing the QS space-time correlation function in Eq. 
( fT5] l with the QS momentum correlation function in Eq. (|4]l, 
one may see a peculiar symmetry: one is transformed to the 
other by the interchange of the emitters and detectors i^. 
Thus the space-time correlations yield the momentum picture 
of the source while the momentum correlations provide the 
information about the source space-time characteristics. 

Some historical remarks are appropriate here. The anal- 
ogy of QS momentum correlations with the HBT space-time 
correlations was first mentioned in paper I2j]. However, not 
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stressing the differences, this^aper triggered a number of mis- 
leading statements, such as [6]: "The interest to correlations 
of identical quanta is due to the fact that their magnitude is 
connected with the space and time structure of the source of 
quanta. This idea originates from radio astronomy and is the 
basis of Hanbury-Brown and Twiss method of the measure- 
ment of star radii". To clarify the situation, Kopylov and 
Podgoretsky wrote a special paper in which they clearly 
stressed the difference between the momentum and space- 
time QS correlation measurements. Particularly, they pointed 
out: "when any of the time parameters characterizing radi- 
ating system becomes very large, the possibility to measure 
the system dimensions practically vanishes since the interfer- 
ence effect remains only in the unobservable small region of 
the energy difference qQ. On the other hand, in astronomy, 
it appears to be possible to measure angular dimensions of 
stars despite their lifetimes can be considered infinitely large". 
Unfortunately, this clarifying paper missed the attention of 
a number of experts in the field of interference correlations. 
Thus even in reviews on the subject one can meet the incor- 
rect statement that there is no principle difference between 
QS correlations in particle physics and astronomy, that both 
are momentum correlations allowing one to determine spatial 
dimensions of the emitting object. An example of such incor- 
rect view of QS correlations in astronomy is given in figure 1 
ifssll . Another example of this widespread error is chapter 1.1 
in a review 15411 . 




star 





correla- 




tion 



FIG. 1 : An example of the incorrect view of spectroscopic QS cor- 
relations in astronomy ^5^1 . Neither the extremely small three- 
momentum difference q between the photons from the same emitter 
(a negligible part of the photon pairs) nor the one between the pho- 
tons from different emitters can directly be measured in astronomy. 
In fact, this figure would be a correct view of femtoscopic QS corre- 
lations if a distant star were substituted by a nearby object emitting 
identical particles at a characteristic space-time distance of the or- 
der of femtometer. 



IV. FEMTOSCOPIC FSI CORRELATIONS 

It can be shown fl3, [H, HI] that the effect of FSI man- 
ifests itself in the fact that the role of functional basis, 
which the asymptotic two-particle state is projected on, is 
transferred from plane waves exp(— //Jijci — ip2X2) to the 

Bethe-Salpeter amplitudes *I'p[p2(xi,X2) = 'Pplpj (xi ,^2) = 
exp(— /PX)\(/~^'*(Ax), where Ajc = jci — X2 = {t,r} is the rela- 
tive four-coordinate, q = q — P{qP)/P^ is the generalized rel- 
ative four-momentum, P = p\+ pi and qP — mi^ — m2^. 

To simplify the calculation of the FSI effect, the Bethe- 
Salpeter amplitude describing two particles emitted at space- 
time points Xi = {f,,r,} and detected with four-momenta pi 
is usually calculated at equal emission times in the pair cm. 
system; i.e. the reduced non-symmetrized Bethe-Salpeter am- 
plitude \\i^^\Ax) is substituted in the two-particle cm. sys- 
tem, where P = 0, ^ = {0,2k*} and Ax = {f*,r*}), by a sta- 
tionary solution V(/L^l (r* ) of the scattering problem having at 
large distances r* the asymptotic form of a superposition of 
the plane and outgoing spherical waves (the minus sign of the 
vector k* corresponds to the reverse in time direction of the 
emission process). This equal time approximation is valid on 
conditions |15, 55t] |f*| <C '«2.i'"*^ for sign(f*) — ±1 respec- 
tively. These conditions are usually satisfied for heavy par- 
ticles like kaons or nucleons. But even for pions, the f* = 
approximation merely leads to a slight overestimation (typi- 
cally < 5%) of the strong FSI effect (see figure 2 and pS!]) 
and, it doesn't influence the leading zero-distance (r* <C |a|) 
effect of the Coulomb FSI. 

Note that in the case of small k*, we are interested in, 
the short-range interaction is dominated by central forces and 
s-waves so that, neglecting a weak spin dependence of the 
Coulomb interaction, the spin dependence of the two-particle 
amplitude enters only through the total spin S. 

On the assumptions (i)-(iv), the two-particle correlation 
function then reduces to the square of the two-particle wave 
function averaged over the distance r* of the emitters in the 
two-particle cm. system and the total spin S of the pair: 



^{puP2) = {\V^^Ar*)\' 



(17) 



For identical particles, the amplitude in Eq. ( fTTl l enters in a 
symmetrized form: 



(18) 



The two-particle approximation in (i) and FSI factoriza- 
tion in the Bethe-Salpeter amplitudes of the elastic transitions 
1 + 2 ^ 1+2 implies a long FSI time as compared with the 
characteristic production time, i.e. the channel momentum k* 
much less than typical production momentum transfer of hun- 
dreds MeV/c. In fact, the long-time FSI can be also factor- 
ized in the inelastic transitions, 1+2^3 + 4, characterized 
by a slow relative motion in both entrance and exit channels 
lf55il56ll . The necessary condition is an approximate equality 
of the sums of particle masses in the channels 1+2 and 3+4. 
In the presence of such transitions the two-channel scattering 
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FIG. 2: The FSI contribution to the 3l"3t*' correlation function calcu- 
lated for the pair velocity v = 0.3,0.5,0.7,0.9 c in a model of inde- 
pendent single— particle emitters distributed according to a Gaussian 
law with the spatial and time width parameters rg = 2 fm and Xq = 2 
fra/c. The exact results (solid curves) are compared with those ob- 
tained in the equal-time approximation (dash curves). 



problem has to be solved and both elastic and inelastic transi- 
tion amplitudes should be taken into account in the averaging 
in Eq. ( fTTb . In practice, the particles 1 , 3 and 2,4 are members 
of the same isomultiplets (as, e.g., in the transition %^ p K^n 
or K^K — > K^K^) so that one can assume the same weights 
and same r* -distributions for the channels 1+2 and 3+4. 



In heavy ion collisions, the effective radius ro of the emis- 
sion region can be considered much larger than the range of 
the strong interaction potential. The short range FSI contribu- 
tion to the correlation function is then independent of the ac- 
tual potential form ifTsilSTll . At small Q — 2k*, it is determined 
by the s-wave scattering amplitudes f^{k*) at a given total 
spin S scaled by the radius ro |15J. For two-nucleon systems, 
the scattering lengths f^{0) are large (up to ^ 20 fm) and this 
contribution often dominates over the effect of QS. For two- 
meson or meson-baryon systems, the scattering amplitudes 
are usually quite small (< 0.2 fm) and the short range FSI 
contribution (including the contribution of the coupled chan- 
nel which is quadratic in the amplitude of the corresponding 
inelastic transition) can be often neglected. This contribution 



cannot be however neglected for the ^T^-system due to rather 
large s-wave KK scattering length dominated by the imagi- 
nary part of ~1 fm gene rated by near-threshold resonances 
/o(980) and ao(980) lH]. It has been recently shown that the 
neglect of the FSI contribution in the analysis of the two-K^ 
correlation function in AuH-Au collisions at V%Af — 200 GeV 



would lead to a noticeable (^ 
relation radius IdSII . 



25%) overestimation of the cor- 



V. FEMTOSCOPIC CORRELATIONS IN BETA-DECAY 

Let us now consider beta-decay of a nucleus with charge 
number Zq, four-momentum po, helicity A-o to a nucleus with 
charge number Z, four-momentum p, helicity X, an electron 
(positron) with four- momentum p^, helicity A-e and an antineu- 
trino (neutrino) with four-momentum pv, helicity A-y. Taking 
into account the point-like character of the weak interaction, 
the equal emission times of the decay particles and the fact 
that the cm. system of the electron and final nucleus practi- 
cally coincides with the rest frame of the initial nucleus (i.e. 
xz = xzo = and k* = pe), one can write the differential decay 
rate in the form: 

• / d^x 1 (x; A's) exp(/pvx)\(/*i^l* (x) 
= Lvs / (^^k*(i^pv5((Bo - CO - cOe - cov) / d^xd^x'- 
■r (x; X's)r* (x'; A's) exp[/pv (x - x')] vL+k** (x)^^^' (x'), (19) 

where the amplitude T(x;A's) is basically determined by the 
distribution of the decaying neutron (proton) within the nu- 
cleus. 

It is instructive to consider the hypothetical situation when 
the energy release in the decay is large and additional particles 
are emitted. Then one could neglect the energy-momentum 
conservation and get, after the integration over the neutrino 
three-momentum (leading to delta-function 5^ (x — x')), simi- 
lar result as in the case of multiparticle production on condi- 
tions (i)-(iv): 



d^k* 



fd\J^ T(x;A's)vL+k':(x) 



¥^;Ux) 



(20) 

The actual energy release in beta-decay is however 
very small so that the integration over the neutrino three- 
momentum does not lead to the diagonalization of the spa- 
tial density matrix ^;)j^/s'r(x;A's)'r*(x';A's). In fact, in so- 
called allowed decays, the neutrino plane wave exp(/p vX J can 
be even substituted by unity. Nevertheless, similar result as in 
Eq. ( l20l i has been obtained by Fermi [34] due to the fact that 
the relativistic Coulomb wave function v'-k** (x) of the elec- 
tron (positron)-nucleus system changes very little within the 
nucleus and can be taken out of the integral in a form of the 
so-called Fermi function F{k* ,Z,R). Neglecting the neutrino 
mass and the nucleus recoil energy (i.e. putting cOy = |pv | and 
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1 



(0 ~ M), one can write in the nucleus rest frame (cOo — Mq): 

,3 ^2 
= 4KF{k* ,Z,R){Mo - M ~ (x),fY, / d\T{x;X's) 



where 



(21) 



F(r,Z,/;) = (|vW(r*)|2)-|vWWl' 



r[ = ^Ze^&e/k* for electrons (positrons), o = (1 — Ze^)'/^ — 
1 . The substitution of the separation r* by the nucleus radius 
R in the second equality in Eq. ( |22] | is justified due to a weak 
r* -dependence of the wave function within the nucleus ; 
the third equality neglects the screening effect of the atomic 
electrons. 

One may see that the nucleus radius enters in the Fermi 
function through the factor (2k*R)^° which is essentially dif- 
ferent from unity only for sufficiently large charge numbers Z. 
At small Z-values, o = and the Fermi function reduces to the 
Coulomb penetration (Gamow) factor Ac (rj) = |v|/*Jkl(0)p = 
27tr|[exp(27tr|) — 1]^^ The sensitivity of the Fermi function to 
the nucleus radius is demonstrated in figure 3. 



VI. CONCLUSIONS 



We have considered femtoscopic QS and FSI momentum 
correlations in multiparticle production and beta-decay, as 
well as spectroscopic QS space-time correlations in the de- 
tected radiation from a distant source. We have demonstrated 
the orthogonal character of the femtoscopic and spectroscopic 
correlations, earlier pointed out by Kopylov and Podgoretsky 
I18J. We have shown that the same functional form of the two- 
particle correlation function in multiparticle production and 
Fermi function in beta-decay (both being equal to the average 
square of the two-particle wave function) is due to different 
reasons. In former case, this result is valid in the approxima- 
tion of independent classical quasi-point-like particle emitters 
and sufficiently small freeze-out phase space density, while in 
the latter case it follows from a weak variation of the elec- 
tron (positron)-nucleus wave function within the nucleus vol- 
ume and a point-like character of beta-decay. It should be 
stressed that a small space-time extent of the emitters alone 
does not guarantee the validity of the approximation of classi- 
cal emitters (the diagonalization of the space-time density ma- 
trix). This approximation may naturally be justified in high- 
energy multiparticle processes due to the minor importance of 
conservation laws in this case. 



Fik*, Z, R) 
25 f 




MeV 



FIG. 3: The Fermi function F{k* ,Z,R) for beta-decay to a final 
nucleus of charge Z = 83 flj a function of the electron ( decreasing 
curves) or positron momentum k* and the nucleus radius R = 2,4,S 
fm ( in decreasing order). 
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